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We investigate the quantum dynamics of Raman-coupled Bose-Einstein condensates driven by 
laser beams that carry orbital angular momentum. By adiabatically eliminating the excited atomic 
state we obtain an effective two-state Hamiltonian for the coupled condensates, and quantization of 
the matter-wave fields results in collapse and revivals in the quantum dynamics. We show that the 
revival period depends on whether the initial nonrotating condensate displays broken U(l) symmetry 
or not, and that the difference may be detected by measuring the motion of quantized vortices that 
are nested in the density profile of the Raman-coupled condensates. We further study the steady- 
state population transfer using a linear sweep of the two-photon detuning, by which the atomic 
population is coherently transferred from an initial nonrotating state to the final vortex state. 
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I. INTRODUCTION 

Controlling and probing of the quantum state of cold 
atoms or molecules are opening up new frontiers for ex- 
ploring the fundamentals of quantum mechanics. At the 
ultracold temperatures characteristic of Bose-Einstein 
condensates (BECs) matter reveals its wave character, 
and the quantum statistics of the atoms or molecules 
play an important role. By controlling both the atomic 
collisions and the properties of laser fields applied to a 
BEC, one may generate tailor made quantum states of 
matter. For example, photons and atoms can interchange 
mechanical momentum, and atoms can therefore be pre- 
pared into a specific center-of-mass momentum state by 
the appropriate choice of spatial light field profile p], Q . 
The light field may also be used for the detection of the 
atomic quantum state 0|. 

A commonly applied notion for atomic BECs is that 
of Bose broken symmetry (BBS) [3, H, El @L or broken 
U(l) symmetry, whereby the state vector for the BEC 
is described via a wave packet of states of varying atom 
number N. In this scheme the fact that the atom number 
is not fixed is traded off for the fact that the macroscopic 
wave function for the BEC, obtained as the expectation 
value of the field annihilation operator, acquires a defi- 
nite phase, hence the name broken U(l) symmetry. The 
BBS approach is typified by the use of a coherent-state 
description for the state of a BEC. In contrast, a number 
conserving (NC) approach is also possible, but in that 
case the phase of the corresponding macroscopic wave 
function has no physical significance, and U(l) symme- 
try is preserved. In the absence of atomic collisions, it 
is known that whether or not the U(l) symmetry is bro- 
ken does not affect any physical observables, and whether 
one describes the state of the BEC as a coherent state 
or a number state is largely a matter of calculational 
convenience 0, 0|. Furthermore, for the case involving 
the interference between two BECs, even for an initial 
state with a fixed number of atoms, quantum coherence 



is built up by the measurement process which leads to 
uncertainty in the atom number, and the resulting in- 
terference pattern is indistinguishable from that which 
would have resulted from a coherent state description 

HI. 

However, this equivalence between the BBS and NC 
approaches holds only when atomic collisions can be ne- 
glected. If one takes atomic collisions into account dif- 
ferences between the two approaches appear in the col- 
lapse and revival times that can occur in the interference 
between BECs with different spatial modes, for exam- 
ple in the interference between condensates [TO, EU, a 
double- well system [12j , or in Ramsey fringe- type experi- 
ments [l3|| . Intuitively, the revival period difference arises 
from the fact that the allowed number difference AN be- 
tween the two BECs is different for the two cases: In the 
NC approach AN is always an even integer, whereas in 
the coherent state approach based on BBS A7V can be an 
arbitrary integer. Since atomic collisions lead to matter 
wave phase shifts proportional to the atom number, the 
interference between two BECs, which depends on AN, 
is markedly different depending on whether U(l) symme- 
try is broken or not, and leads to a factor of two difference 
in the revival period between the two approaches. This 
factor of two difference has recently been observed for 
the first time in an experiment that looks at the interfer- 
ence between atoms released from a lattice of double well 
systems [14J: Depending on the preparation conditions, 
the atoms trapped in each individual double well in the 
lattice may be best described as a number state or a co- 
herent state, and this leads to a factor of two difference 
in the revival time in the interference pattern that arises 
when the atoms are released and allowed to interfere. 

Collapse and revivals have previously been investigated 
experimentally in the Jaynes-Cummings model of a two- 
level atom interacting with a quantized laser field [l5[, 
wave packet dynamics of Rydberg atoms [l6[, matter 
waves in an optical lattice [17[, and most recently for 
a matter wave in a lattice of double wells [l4[. Never- 
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theless, the use of collapse and revivals to test the fun- 
damental notion of U(l) broken- symmetry as applied to 
a single BEC via the period of the collapse and revivals 
has not been realized experimentally so far. Walls and 
co-workers [10|, showed that this could be realized by 
coupling an initial BEC with a given mode structure to 
a second BEC mode, thereby realizing a two-mode sys- 
tem. The period of the revivals of the system would then 
act as a test of whether or not the initial isolated BEC 
exhibited BBS or not. Search |13| has proposed precisely 
such an experiment based on a Ramsey fringe approach, 
but this has not been realized to date. 



In this paper, we consider the case of Raman-coupled 
BECs using Laguerre-Gaussian (LG) beams [18| as a po- 
tential testing ground for the BBS description of an initial 
BEC. Because of the helical phase structure of the laser 
fields, orbital angular momentum (OAM) is transferred 
from the light to the atoms and results in a BEC in a 
coherent superposition of two components with distinct 
center-of-mass OAM 0,[l9|,[2O|. (In another scheme for 
generating vortices a de-centered Gaussian beam is ro- 
tated around the center of a trapped BEC [2l[, but we 
do not analyze that case here.) By considering BECs 
with relatively small numbers of atoms, and quantizing 
the matter- wave fields, the system can be mapped onto a 
Hamiltonian that describes the dynamics of an ideal two- 
mode system. Then the granular nature of the matter- 
wave field becomes important, and we explore the col- 
lapse and revivals in the system both with and without 
broken U(l) symmetry. In particular, we show that the 
difference between the two approaches is directly visible 
in the motion of quantized vortices that appear in the 
density profile of the Raman-coupled condensates, the 
density profile being a relatively easy physical quantity 
to observe. To put our proposal in context, collapse and 
revival dynamics have previously been observed using the 
spatial density profile in a two-component 87 Rb BEC [22[, 
and the Raman vortex coupler using LG beams has now 
been realized experimentally [23], both of these experi- 
ments being performed for large atom numbers where the 
granular nature of the matter wave fields is not relevant. 



This paper is organized as follows. In Sec. [II] we 
formulate the problem, and map the three-level A-type 
atomic system coupled with LG fields onto a two-mode 
Hamiltonian in the angular-momentum representation. 
In Sec. HIH we focus on the collapse and revival dynam- 
ics for initial number states and coherent states for the 
BEC, and show that these may be measured using the 
spatial density profile of the Raman-coupled condensates 
as well as the atom statistics. Section [TV] summarizes 
the results in this paper, mentioning the time-dependent 
two-photon detuning for the stimulated Raman popula- 
tion transfer. 



II. FORMULATION OF THE PROBLEM 

In this section we formulate the two-mode approxima- 
tion for two Raman-coupled BECs and introduce the an- 
gular momentum representation used to solve the quan- 
tum dynamical system. 



A. Hamiltonian 

We consider a quantum-degenerate sample of ultracold 
bosonic atoms with three participating levels |a), |6), and 
\e) arranged in a A configuration, see Fig.[TJ The sample 
is irradiated by a pair of laser fields of frequencies cji 5 2, 
transitions between \a) and \b) being dipole forbidden. 
The system Hamiltonian H = H\ + Haf is given by [24| 
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where are creation and annihilation field opera- 

tors for bosonic atoms in the states j = a, 6, e. Atomic 
collisions are included via the pseudopotential coefficients 
Xjj' = 4tt^ 2 ajj'/M, where a,jj> is the s-wave scattering 
length between two atoms in the states j and j', and M 
is the atomic mass. The coupling of the atoms to the 
classical optical fields is described by the Rabi frequen- 
cies •> wnere i denotes the positive and negative fre- 
quency components. The atoms are further supposed to 
be confined in harmonic traps, so that the center-of-mass 
Hamiltonian for atoms in each state j = a, 6, e is given 
by 
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the oscillator frequencies along the 



radial (r) and longitudinal (z) directions for the different 
atomic states. 

If the detuning A from the excited state \e) is suffi- 
ciently large (see Fig. 1), the population of that state 
remains negligible and spontaneous emission may be ne- 
glected. The excited state can then be adiabatically elim- 
inated as 
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FIG. 1: A atomic level scheme: Three hyperfine states \a), \b), 
and \e) are coupled with light fields of Rabi frequencies Qi, 
Q2, and laser frequencies cji, 002, respectively. The common 
detuning A is sufficiently large that the excited state |e) is 
not significantly populated. 



and the Hamiltonian for the atom-field interaction re- 
duces to 



-2ft 



/ 



cPr 



l# ) fti" ) f i(aJ2 - Wl) ** t 



h.c. 



,(6) 



where we have omitted the diagonal terms due to the 
ac-stark shift, which create an effective potential for the 
atoms but do not affect the dynamics. The atom- field 
interaction Hamiltonian Haf in Eq. ([6]) reflects the fact 
that if an atom makes a transition from state \a) to the 
excited state via absorption of a photon of frequency uj\ , 
this will be rapidly followed by emission of a photon of 
frequency 002 and a Raman transition to the state \b). 

Bringing these results together, the system Hamilto- 
nian H = Ha + Haf is now given by 



H A = 



Haf = 
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This Hamiltonian is the starting point for our analysis. 



B. Orbital angular momentum of light 

Light has two kinds of angular momentum: a spin- 
angular momentum associated with its polarization, and 
an orbital angular momentum (OAM) associated with 
mechanical rotation [25[. The system described by the 
Hamiltonian ([7j) corresponds to a familiar coherent cou- 
pler between the atomic BECs in the states \a) and \b) 
if both light fields are plane waves or Gaussian beams, 
that is, beams that carry no OAM. Then the absorption- 
emission process transfers the atom from the state \a) to 
the state \b) that is shifted in momentum by the two- 
photon recoil. On the other hand, if one or both of those 
laser fields carries OAM, the light field can impart its 
OAM to the atoms, and hence angular- momentum trans- 
fer occurs associated with the coherent coupling [26[. In 



this case, the two-photon recoil is associated with the 
difference of angular momenta of two light beams, and 
one realizes a Raman vortex coupler. The Raman vortex 
coupler has previously been studied theoretically in the 
limit of large atom numbers where a treatment based on 
coupled Gross-Pitaevskii equations is applicable [18], and 
an experimental realization of this case has been reported 
[23( | . Here we explore the opposite limit of low atom num- 
bers, where the granular nature of the quantum matter 
wave fields becomes relevant. 

A typical set of light fields that carry OAM are 
Laguerre-Gaussian beams whose mode functions at the 
beam waist z = are given by [25[ 
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where wq is the focused spot size at z = 0, p is the radial 
mode number, I the winding number, which describes the 
helical structure of the wave front, and 
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are Laguerre polynomials. By using LG modes for 
the Raman coupling between the matter wave fields in 
the Hamiltonian ([7j) we can therefore create a matter- 
wave vortex even though the initial state is nonrotating. 
Henceforth, we consider only the LG modes with p = 0, 
and we thus need only consider one quantum number Z, 
to specify the OAM. 



C. Two- mode Hamiltonian 

We assume that all atoms are initially in the state \a) 
with total angular momentum zero, and that this state is 
Raman-coupled to the state \b) with a winding number 
I by using a pair of LG laser fields, one with a wind- 
ing number zero and the other with winding number I. 
We further assume that two-body interactions are weak 
enough, and the external harmonic trapping frequencies 
w r j = oo r and u z j = uo z are large enough, that the mode 
profiles of the two atomic states are not significantly mod- 
ified from their single particle forms i/)o,i(r) [12|, [27[. 

Under these assumptions the matter- wave field opera- 
tor may be expressed as 



$(r,t)=a(^(r)+i(#W, 



(10) 



where a and b are the bosonic annihilation operators for 
the states \a) and |6), and the positive frequency compo- 
nents of the Rabi frequencies are 

Sl[ + \r) = n 1( f>i(r) , n ( +\r) = n 2 <f>o(r). (11) 

When the LG beam remains well collimated over the lon- 
gitudinal extent of the condensate, d z <C kow 2 /2, the 
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laser mode functions 0o,z at the beam waist and the con- 
densate mode functions ^o,z may be expressed as 

<j>l{r) = -L ( — ) e -r 2 /2wl e il^ e ikz^ (12) 

*<" = 7lR(i) ,, ^ e °'** w - (13> 

We assume that Z > 0, and that the radial and lon- 
gitudinal atomic oscillator ground state widt hs are the 
same for both atomic states d r = ^/h/(Muj r )^ and d z = 
yjhj (Mco z ) [Hf. We further assume that the longitudi- 
nal trapping is much stronger than the radial trapping, 
d r > d Zl and that the longitudinal wave function remains 
fixed as the ground state ip z (z) = e~ z ' l / 2d z / ' ^/^J^d~ z . 
The two- mode approximation is then valid under the con- 
ditions Ndij < dk [12], where i,j — a,b and k — r,z. 

Substituting Eqs. ([T0|) - (fT2]) into the system Hamil- 
tonian Q with Ej — j d 3 ripj Ho^j: and transforming 
to a rotating frame using the unitary transformation 
U = exp{-i[E a (tfa + ttb)/H - (uo 2 - u)i)tt b]t}, finally 
yields the reduced two-mode Hamiltonian 



H — fajJ r 



(l + 5/u r )Pb + G aa a) a) da + G bb b r b r bb 
+ 2G ab d r b r ba - g(a r b + atf)] , (14) 



where we have defined S — 002 — oo\ + (E b — E a )/H as 
the two-photon detuning, which may be controlled ex- 
perimentally. The dimensionless coupling constants for 
the atom-atom and atom-light interactions are given by 
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(16) 



D. Angular-momentum representation 

In order to study the quantum dynamics described by 
the Hamiltonian ([T4]h we introduce the following opera- 
tor combinations that map our problem to Schwinger's 
angular- momentum representation [l2[ , 
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Jz 



db\ J_ = 
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These operators obey the SU(2) commutation relations 
[Ji, Jj] = itijkJk where e^k denotes the Levi-Civita an- 
tisymmetric symbol, and the Casimir invariant is found 
to be J 2 = J 2 + Jy + J z = j (j + 1) with j = N/2. For 
the basis |j, m) with m = — j, — j + 1, . . . , j, which are 
eigenstates of J Z1 the operators satisfy 
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J 2 \h m) 
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j(j + l)|j,m), (24) 
y/(j±m+ l)(j =F m)\j, m). (25) 



In the angular-momentum representation the Hamilto- 
nian ([T4]) then becomes 



H = h z J z + X J z -2gJ x 



where 
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The angular basis states | j, m) and the Fock basis states 
l^a, ^6)f are related by \j, m) = \N/2 — m, N/2 + m)p. 



with 
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III. QUANTUM DYNAMICS 

In this section we explore the quantum dynamics of 
the two-mode model, first with no atomic collisions, and 
then we study the collapse and revivals that occur with 
atomic collisions. In particular, we highlight that the 
collapse and revivals may be monitored by following the 
motion of quantized vortices that appear in the spatial 
density profile, and we shall also study the resulting atom 
statistics. 



The parameter £ = d r /wo characterizes the ratio of 
the radial oscillator ground state width to the focused 
light spot size, and the Kronecker delta enforces angular- 
momentum conservation for the atom-field interactions. 
Henceforth the energies, lengths, and angular momenta 
are measured in units of hw ri d r , and H, respectively. 



A. No atomic collisions 

In the absence of atomic collisions, Gij = 0, giving 
h z = 1 + 5 and x = 0, the two- mode Hamiltonian reduces 
to the form of two simple coupled harmonic oscillators. 
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Then the solution of the Heisenberg equations of motions 
for a(t) and b{t) is found to be 
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and the effective Rabi frequency v is defined as 
_ N /(2 5 ) 2 + (I + 5? 



(29) 
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We consider an initial state where all atoms are con- 
densed in the hyperfine state \a) with zero OAM. If we 
take a number conserving initial state with definite atom 
number, \V(0)) N = \N, 0) F = |7V/2, -N/2), then the half 
population difference J z = (b^b — a) a)/ 2 evolves in time 
as 



n(*(0)\JzW))n 



(Z + <y) 2 + (2g) 2 cos(2i/t) 
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The total-angular-momentum operator is just propor- 
tional to the number operator of the state \b), 



L = lb r b = l(J + J z ), 



(34) 



and the fluctuations in angular momentum coincide with 
the number fluctuations of the state \b). The time evo- 
lution of the angular momentum for an initial Fock state 
is given by 
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We note that the maximum population transfer occurs 
at the resonant two-photon detuning 5 = — /, and the pa- 
rameter g in Eq. JTBJ) describing the single-particle Ra- 
man coupling between the LG field and the BECs de- 
creases with increasing winding number I. 

We now repeat the same calculation, but using a 
Glauber coherent state as an example of an initial BEC 
displaying BBS 
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with \a\ 2 = N being the average number of atoms. 
The number and phase fluctuations of the coherent state 



are given by An = |a|, and ~ 1/(2 |a|), respec- 
tively. The time evolutions of the half population differ- 
ence a(fy(0)\J z \ty(0)) a and the total angular momentum 
a (^/ (0)\L\^/ (0)) a are easily shown to coincide with Eqs. 
([33]) and (|35|) with the replacement N — > N. In other 
words, in the absence of the atomic collisions, the phys- 
ical quantities are therefore independent of whether the 
state of the initial BEC exhibits BBS or not. 



B. Collapse and revivals 

The two-mode Hamiltonian exhibits a rich variety of 
quantum dynamics depending on the relative magnitude 
of the atomic collisions and atom-field coupling, and the 
OAM of the LG beam. The system dynamics can be 
categorized into three typical parameter regimes [29^ 
namely (i) the Rabi regime g ^> \N , (h) the Joseph- 
son regime g <C ^ 9^ 2 i an d (hi) the Fock regime 
gN <C X- Oscillations between the two BEC states are 
almost perfect in the Rabi regime at the resonance condi- 
tion h z = 0, but the transfer between the coupled BECs 
reduces as the Josephson regime is approached. In the 
Fock regime, the Rabi oscillations are suppressed and the 
number in each internal state is almost a constant, like 
the case of self-trapping in a double- well potential or the 
Mott state in optical lattices. The critical point of the 
self-trapping transition is xN — 4g, which lies in the 
Josephson regime (ii) . Here we concentrate on the situa- 
tion where the system is in the Rabi regime, considering 
for concreteness the parameter values N = N = 100, 
X = 5 x 10 -3 , and g = 10, corresponding to a conden- 
sate with a mesoscopic total number of atoms. In the 
following numerical study we choose I = 2 as a repre- 
sentative example, and take the two-photon detuning as 
6 = -1 + (N - l)x-, so that h z = 0. With these pa- 
rameters we anticipate almost complete Rabi oscillations 
of the population difference between the two condensate 
components. 

To proceed we introduce two initial states that repre- 
sent NC and BBS states for the initial BEC. First, for 
the NC state we choose a state with a well-defined atom 
number \N, 0)f in Fock-space representation. The subse- 
quent quantum dynamics of the coupled BECs may then 
be represented in the angular-momentum basis as 
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^2 A ™(t)\j,m), 
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where j = iV/2, and the normalization is given by 
l^m(^)| 2 = 1- On the other hand, for an initial 
condition with BBS where the initial state is a coher- 
ent state \a,0), the time evolution of the coupled BECs 
may be represented in the angular-momentum basis as 
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FIG. 2: Time evolutions of (a) (J*)jv, (b) (J z ) a , (c) (J*)iv, 
(d) (J x ) a for AT = JV = 100, x = 5 x 10" 3 , and g = 10. 



period is also found in the double- well system and for the 
Ramsey fringe experiment proposed by Search, and the 
revival period agrees with that he found for the resonant 
case [13| . 

For the sake of completeness we comment briefly on 
the effect on our findings of selecting a nonzero value of 
the parameter h z . For the case h z ^ 0, which would be 
possible using a Feshbach resonance, the system enters 
the Josephson regime as \- increases, and this reduces 
the collapse time for both the Fock and coherent states. 
The Rabi oscillations are then no longer perfect, and the 
coherent state collapses much faster than the Fock state. 
Although the revival time becomes longer, the factor of 
two difference between Fock and coherent initial states 
remains. This is different from the situation in the Ram- 
sey technique, where the system evolves freely in time 
(without coupling between the BECs) between the Ram- 
sey pulses and the coupling term proportional to J x in 
the Hamiltonian is zero. In the present case, by contrast, 
the Raman coupling is on at all times. 

Finally we emphasize that the factor of two difference 
in revival period depends only on whether the initial state 
is NC or not. For instance, if the initial state is chosen as 
number-conserving coherent state with binomial statis- 
tics |3l|, then the period of the revival will still be half 
that of the Glauber coherent state. 



where there is no dynamical coupling between subspaces 
with different total number of atoms since the Hamilto- 
nian conserves the atom number [301 ]. 

Figure p] shows the numerically calculated time evolu- 
tion of (J x )n and (J z )n for an initial Fock state, and of 
(Jx)a an d (J z )a fc> r an initial coherent state. In both 
cases atomic collisions cause the anticipated complete 
Rabi oscillations between the two BECs to be modulated 
by a sequence of collapses and revivals, the revivals be- 
ing a characteristic consequence of the granular nature of 
the matter wave fields. The components J XjV character- 
ize the coherence between the states \a) and |6), while J z 
corresponds to half their population difference. We note 
that (J x ) exhibits a revival in the middle of the collapse 
of (J z ) for both coherent and Fock initial states. 

There is a significant difference between the period of 
revivals for the NC and BBS initial states. In the Rabi 
regime, and for the resonant case h z = 0, the collapse 
time and revival time in the population difference (J z ) 
are numerically found to be 



C 



collapse xv /]v' 

for the initial Fock state, and 
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for the initial coherent state, respectively, where the con- 
stant is C ~ 10. This factor of 2 difference in the revival 



C. Quantized vortex motion 

The results of the previous section show that the col- 
lapse and revivals in the dynamics of the quantum vortex 
coupler can be used to test whether the initial state of the 
BEC is described as a NC state or a BBS state, and these 
results are in perfect accord with previous works on this 
topic. However, measuring the components of the angu- 
lar momentum in Schwinger's representation raises issues 
of how to detect these quantities. Here we propose to de- 
tect the collapses and revivals by monitoring the spatial 
density profile of the Raman-coupled BECs, and in par- 
ticular, the motion of the quantized vortices that appear 
in the profile. Collapse and revival dynamics have pre- 
viously been observed by monitoring the spatial density 
profile in a two-component 87 Rb BEC [22j], and a Raman 
vortex coupler using LG beams has now been realized 
experimentally (23[. These experiments were performed 
for large atom numbers where the granular nature of the 
matter wave fields is not relevant, but they show that the 
proposed approach has some validity. 

To proceed we investigate how the quantum dynamics 
can manifest itself in the spatial density profile of the cou- 
pled BECs. From Eq. ([TO)) , the density operator p(r,t) 
may be written as 

p(r, t) = t(r, t)*(r, t) = p (r, t) + pi(r, t), (42) 




FIG. 3: (a) Time sequence of density profiles (p(r,t))jv given by Eq. ([42]) . and (b) expectation values of the pseudospin 
components during one Rabi cycle. Solid curve with the filled dots shows {J z )n, dotted curve with the open dots {J v )n, and 
dotted curve with the filled dots {J x )n, respectively, (c) Off-diagonal part of the density (pj(r,t))jv given by Eq. ([44]) . Upper 
and lower panels show the typical patterns for the first half and latter half of a Rabi cycle, respectively. The alignment of the 
two density maxima (mimina) rotates by tt/2 anytime when the {J v )n crosses the zero. 



where 
Po(r,t) 



\Mr)\ 2 

r l \Mr)\ 2 r 
r l \Mr)\ 2 



l + ^)j(t)-(l-^)j s {t) 



(43) 

j_(t)e ilc ? + j+(t)e- il( P 
J x (t) cos(l(p) + J y (t) sm(l(p) , (44) 



and the total density profile for both atomic states is 
given by (p(r,t)) = (po(r,t)) + (pi(r,t)). 

From this equation we see that the spatial density 
profile is directly related to the pseudospin components, 
hence it will also display the collapse and revival phe- 
nomenon. Since the qualitative properties of the density 
profile are the same for both the Fock and coherent initial 
states, with the caveat that the revival period is twice as 
long for the coherent state, here we show results for the 
case of an initial number state. 

Figure 02(a) shows the evolution of the density pro- 
file (p(r,t))N over one Rabi oscillation for times suf- 
ficiently short that the first collapse is barely notice- 
able, and Fig. (3jb) shows the corresponding evolution of 
(Ji)N/N,i = x,y,z for the number state here times are 
scaled to the revival time T r ^ val . For the earliest time 

V^revL = 5 X 10 " 4 in theSe P lot8 > (Jz)n/N ~ -1/2, 

meaning that the condensate is almost entirely in the 
nonrotating Gaussian mode, and the left-most gray-scale 
density profile in Fig. [3](a) shows this Gaussian density 
profile. As time progresses, however, the Raman cou- 
pling transfers more population to the rotating BEC 
mode with I = 2, and the second through fourth den- 
sity profiles in Fig. [3](a) show this process in which 
two vortex cores with winding number I = 1 progres- 
sively make their way towards the center of the con- 
densate. Halfway through the Rabi oscillation the BEC 
is composed purely of the rotating state, corresponding 



to the time t/T^ ival = 6 - 25 x 10 ~ 4 in Fi S- G3( b ) where 
(Jz)n/N ~ 1/2, and the density profile is now that for 
a pure matter-wave vortex with I = 2 centered at the 
origin. For longer times the process reverses up until 

^/ ^revival ~ ^ x 10 -4 when the next Rabi oscillation 
starts. 

During a single Rabi oscillation the distance between 
the two I = 1 matter-wave vortices is related to the pop- 
ulation difference (J z )n- The angle of the line joining the 
two I = 1 vortex cores (where the density goes to zero) 
may be understood by looking at the off-diagonal part 
of the density profile (p/(r, As seen from Eq. (|44]h 

the angle of the line joining the vortex cores is deter- 
mined by the pseudospin components (J x ) and (J y ). For 
example, the spatial distribution of the off-diagonal part 
of the density profile at the angle <p = is always pro- 
portional to (J x ), whereas at the angle tp = 7r/4 it is 
proportional to (J y ). For the parameters of this example 
we have that \(J x )n\ <C \{J v )n\^ see Fig. [3fb). Then, 
viewing the expectation values of the pseudospin compo- 
nents as parameters on a Bloch sphere, one Rabi oscilla- 
tion corresponds to a round trip between the south pole 
({Jz)n/N = -1/2) and the north pole ({J Z ) N /N = 1/2), 
with the shortest path being taken on the Bloch sphere. 
Since the component (J x ) is small compared to the other 
two components, the off-diagonal term is approximately 
given by (pj) ~ r l \ifjo(r)\ 2 (J y (t)) sin (lip)/ y/l\, which illus- 
trates the change in profile in Fig. G2c) for times before 
and after the system evolves halfway through the Rabi 
oscillation in Fig. [31(b). When the value (J v )n crosses 

zero halfway through the Rabi oscillation at t/T^j ^ = 
6.25 x 10 -4 in Fig. [31(b), the positions of the off-diagonal 
density maxima and minima are interchanged. This gives 
rise to the change in the angle of the line joining the two 
I = 1 vortex cores in Fig. [31(a). We remark that if the 
dynamics takes the longer path on the Bloch sphere, i.e., 
as happens when J x and J y become comparable and ex- 
change their values, the position of the vortices can rotate 
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FIG. 4: Density profiles (p(r,t))jv (upper panels) and half 
population difference {J z )n/N (lower panel) for times allow- 
ing for collapse and revivals. The density plots are shown at 
the times indicated by small circles in the lower panel. 



around the condensate. Hence the radius and angle of the 
vortex position generally follow a trajectory on the Bloch 
sphere. 

Figure [H shows a time sequence of density plots over 
time scales long enough to allow for a collapse and re- 
vival, the lower plot showing the evolution of (J z )n/N 
versus time. For short times the system undergoes Rabi 
oscillations as in Fig. [3l but modulated by a collapse en- 
velope are shown for times t/T^ val < 0.1, in the upper 
left panel. In the collapse region 0.1 < t/T^ val < 0.9 
when {J z )n/N — > 0, the density profile is to all intents 
and purposes stationary and the vortex cores assume 
fixed positions, see the central upper panel. Finally, for 

'^revival > O-^ the & TS ^ revival starts, the vortex cores 
get back into motion, and the Rabi oscillations start 
again, see the right most upper panel. Thus the mo- 
tion of the vortex cores in the density profile provides a 
convenient means to monitor collapse and revivals in the 
quantum dynamics of the matter wave vortex coupler. 
Density measurements thus represent a powerful tool to 
measure the dynamics of the collapse and revivals as well 
as the motion of the condensate on the Bloch sphere, and 
as such allow one to measure the factor of two difference 
in the revival times for initial state displaying Bose bro- 
ken symmetry or being number conserving. 
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FIG. 5: Mandel's Q parameter defined by Eq. (|45|) for initial 
(a) Fock and (b) coherent state. The value becomes —1 near 
the revival for the initial Fock state, while it is always larger 
than —1. 



state [Fig. El (a)], Q a = —1 initially, and Q a (t) returns to 
that value periodically near the revivals. For an initial 
coherent state, on the other hand, the Mandel parame- 
ter is initially Q a = 0. It then becomes negative near 
the revivals, indicative of a sub-Poissonian distribution, 
but never reaches the Fock state value Q a = —1. Note 
also that the average values of Q a for the coherent and 
Fock state approach each other as the number of atoms 
increases. 



D. Atom statistics 

For completeness we now briefly discuss the atom 
statistics of the system as it undergoes a series of col- 
lapses and revivals. A familiar measure of the atom- 
number fluctuations is Mandel's Q parameter [32[ 

((An,) 2 ) - (n d ) . 

Qj = 7^ r , J = a,o, (45) 

\ n j) 

where Q = corresponds to a Poissonian, — 1 < Q < 
to a sub-Poissonian, and Q > to a super-Poissonian 
distribution. In particular, Q = — 1 corresponds to a 
Fock state with (An,) = 0. As expected and illustrated 
in Fig. [5l the Mandel parameter also exhibits the factor 
of two difference in the revival times depending on the 
U(l) symmetry of the condensate. For the initial Fock 



E. Effects of thermal dissipations 

To assess the feasibility of an experimental demonstra- 
tion of these predictions, we consider the specific exam- 
ple of a 87 Rb condensate with Raman coupling between 
the two hyperfine states \a) = \F = l,rap = —1) and 
I b) = \F = l : mF = 1). The trapping potentials of these 
two states are almost identical, V\(r) T^(r). It is ex- 
perimentally clearly desirable to have a relatively short 
revival time, so that a tight condensate confinement is 
preferable. Taking u r = 2tt x 100 Hz and uj z — 2ir x 300 
Hz, the corresponding oscillator lengths are then d r ~ 1 
/im, d z c± 0.6 jam. All s-wave scattering lengths are 
approximately the same, a aa c± c± a a & c± 5.5 nm. 
For this set of parameters, the validity of the two- mode 
model requires that the total number of atoms has to be 
N < 200. For the case of an OAM 1 = 1, these values 
result in the dimensionless parameters G aa = 2Gbb = 
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FIG. 6: Population transfer from the state \a) to \b) with the 
linear sweep of the two-photon detuning. The insets and the 
numbers express the intensity of the interference pattern and 
the corresponding (dimensionless) time. 



2G ab ~ 3.7 x 10~ 3 and thus x = G aa /2 = 1.8 x 10~ 3 
[see Eqs. ([T5]) and ([T7|)]. Typical time scales are then 
estimated as T Rabi ~ 1.4 x 10" 3 s, T c ( o ^ apse ~ 0.5 s, 

Revival - 5 - 5 s > and Revival - 1L0 s > respectively. Hence 
the difference in the revival periods of the Fock and co- 
herent states should be observable in a condensate with 
a relatively long lifetime on the order of a few tens of 
second. 

We note that there is a caveat with respect to our 
proposal in that it involves multiply-charged matter- 
wave vortices with winding numbers whose magnitude is 
greater than one. In particular, it is well known that if a 
multiply charged vortex is initiated and allowed to evolve 
freely it will quickly decay into multiple singly charged 
vortices in the presence of even a small background com- 
ponent of thermal atoms [33]. The decay time for a 
doubly-quantized vortex ranges from one to a few tens of 
milliseconds, depending on the total number of atoms. 

However, this does not pose a problem for the Raman 
coupler analyzed here, as the matter- wave vortices are 
not freely moving, but rather are externally driven by 
the LG beams. More specifically, the applied LG laser 
fields continually drive the matter- wave field between the 
zero-charged (Z = 0) and doubly charged (I = 2) vortex 
states with a time period corresponding to the inverse 
Rabi frequency TR a bi ^ 1.4 x 10 -3 s. (This is true even 
in the collapse region where the time-independent density 
profile arises from an average over many oscillating com- 



ponents with slightly different periods.) Thus, as long as 
^Rabi <^ Td there will be negligible decay of the doubly 
charged vortex created in any Rabi cycle before it is re- 
turned to the zero-charged state and the next Rabi cycle 
starts. 

IV. SUMMARY 

In summary, we have studied the quantum dynamics of 
the vortex coupler using the Laguerre- Gaussian beams. 
As a consequence of the quantization of the matter- wave 
field, the atoms are found to undergo a series of collapses 
and revivals whose period is directly observable as an off- 
axis motion of the quantized vortex cores that appear in 
the condensate density. The angle of the vortex also pro- 
vides direct information on its path on the Bloch sphere 
that describes its dynamics in the Schwinger representa- 
tion. An important feature is that the characteristic time 
scale of the collapse and revivals differs by the factor 2 
depending on whether the description of the condensate 
is in terms of a number-conser ving state or a state with 
broken U(l) symmetry 0, El El- 

We finally note that instead of generating the collapse 
and revivals used to test the appropriate theoretical de- 
scription of the condensate, our scheme can also be modi- 
fied slightly to permanently transfer the condensate pop- 
ulation to the state \b). All that is required is to sweep 
the two-photon detuning linearly in time at an appropri- 
ate rate. We illustrate this technique with the parame- 
ters N 100, x = 5 x 10~ 3 , g = 2, and S = 20 - U 
in Fig. [6l Since both the initial and final states are in 
the Josephson regime, no detectable collapse and revival 
occurs in any observables. The J x and J y components 
have maximum values when the detuning 5 crosses the 
resonance, and a significant change in J z occurs at this 
point. Note that a population transfer is achievable inde- 
pendently of the U(l) symmetry, although the variances 
and atom statistics differ significantly in both cases. 
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